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54 WOODWARD. ON THE ERRORS INCIDENT TO VALUES 

ON THE ERRORS INCIDENT TO VALUES IN TERPOLATED BY MEANS OF 

FIRST DIFFERENCES FROM TABLES OF LOGARITHMS, NATURAL 

TRIGONOMETRIC FUNCTIONS, ETC. 

By R. S. Woodward, Wa.shington, D. C. 

The theory of the errors incident to value.s interpolated by means of fir.st 
differences from a table of logarithms, natural sines, etc., appears to be a matter 
of considerable difficulty. Our text-books are notably deficient in information 
on this subject or notably erroneous in such information as they give; indeed 
few writings on the subject are free from error. 

The most common departure from correct procedure in treating this class of 
errors is in assuming the applicability to them of the law of error of least 
squares. In addition, the magnitude and continuity or discontinuity of the 
errors involved are frequently misunderstood. Here also as in other branches of 
the theory of errors, the actual, probable, mean, and average errors are apt to be 
confounded with one another. 

Since attention has recently been called in the Annals by Prof Howe to 
the errors incident to the most frequently employed species of interpolated values, 
it may be of interest to point out some of the essential features of the correct 
theory of such errors. 

In a paper by the writer published in The Analyst, 1882, four species of in- 
terpolated values were considered. Two of these depended on the tabular dif- 
ferences and two on the true differences. It was shown that the latter are less 
precise than the former, or, paradoxical as it may seem, that the use of the true 
differences is less advantageous than the use of the tabular differences. Of the 
other two species, the first retains one or more figures beyond the tabular number 
of figures, and the second (which is most commonly employed) rounds the inter- 
polated value to the nearest unit of the last tabular place. It will be convenient 
to refer to these two as the the first and second species. 

The theory of the errors of the first species has been fully worked out by 
Bessel {Astronomische Nachrichten, Nos. 358 and 359, 1838) and affords a basis 
for the theory of the errors of the second species. 

I am not aware that the theory of the errors of the second class has been 
completely developed. In the paper referred to above, however, I have indicated 
the process of deriving the laws of frequency of the errors of the second species 
and have given the laws and consequent probable errors for several special cases. 
I am now able to extend the process and give more general results. 

The actual error £„ of an interpolated value of the first species is expressed 
by the equation 

£„ = (I -/)£, + /S,. (I) 
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in which s, and s^ ^^^ the actual errors of the tabular logarithms or numbers, 
and / is the interpolating factor, which varies from o to i. e, and % are expressed 
in units of the last tabular place and may have any value from — ^ to + ^ ; 
hence also e„ may have any value between the limits — ^ and + ^. 

The actual error s of an interpolated value of the second species will equal 
that for the first species plus the error which comes from rounding the product 
of the interpolating factor and the tabular difference to the nearest unit. Hence, 
if wc call the error from the latter source £3, we shall have 

£ = (l -/)£, + /£2 +63. (2) 

It must now be observed that £3 differs from Sj and e^ in two important re- 
spects, viz : 1°, it is discontinuous, and 2°, it is dependent on t. Thus for / = ^, 
Sj must be either o or ^. Again, for i = ^,s^ is restricted to one of three values, 
viz: o, + ^ and — ^. On this account the maximum value of e will vary also 
with /. For / = ^the maximum value of e = ^ + ^= i, but for i^=^ the 
maximum of £ = -^ + -J- ^ |^; this maximum cannot always be i as some of the 
books assert. Moreover, if / is small relatively to the largest tabular difference 
for the table used, £, will be correspondingly small. Thus, in a 5-place table of 
logarithms, the largest tabular difference is 44 ; and if / := o.OOi the maximum 
value of £3 is 0.044 and the maximum of £ is 0.544.* 

To determine the probable, mean, and average errors corresponding to £„ and 

£( we must know the laws of frequency of these actual errors. If we represent 

the probability of the occurrence of an error £ by f (s) ds, it may be shown thatf 

for the system of errors represented by equation (i) 

1 I 5 
If (£„) = ?-~-~-Y foi" values of £„ between — \ and — (^ — t), 
(I /) / 

^= ^^ for values of £„ between — (-^ — i) and + (i — i)< (3) 

=: j^ -^ for values of £0 between + (i — ^) and + \. 

As our present purpose is to discu.ss the efrors of the second species, we 
shall not derive the probable and other critical errors of the first species. It is 
sufficient to remark that these critical errors show that an interpolated value of 
the fir.st species is in general much more precise than a tabular value.f 

To derive the laws of frequency of the errors £ of the second species, it is 
necessary to combine the continuous errors £„, whose laws of frequency are 

*'rhese facts would seem to be perfectly plain, but experience justifies the assertion that they are not 
commonly appreciated. 

tSee Bessel's paper or that of the writer referred to above. 
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given by (3), with the discontinuous errors s^ First consider the case in which / 
is a fraction of the form 



n 



in which n is an integer.* Then the possible values of £3 are the following : 
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for n even. 



These errors are all equally likely to occur. The sign of the error \ for n 
even is arbitrary. It will here be considered positive only ; and, in deriving the 
laws of frequency of the errors t for n even, they will be considered with re- 
spect to magnitude simply, or as if they all had the same sign. It will be ob- 
served that the numerical maximum of s is 



2n — I 
2n 



for n odd. and 



-\- \ -^ \= -\r \ for « even. 

Now if we add each of the above n values of % to £„ in (3), we shall have n 
sets of probability functions, <f (£„ + s,) ^ ^ (j). If then the element-areas 
i5 (e) dt defined by these several functions be summed for every value of s, the 
result will be the probability area for the whole .system of errors. The area en- 
closed by the straight lines whose equations are given by (3) and the axis of £(,, 
IS I. Hence the sum of n such unit areas is n, and we must divide the results 
of the above summation by n in order to have the integral of if (s) dt over the 
whole range of error, or the sum of the probabilities of all possible errors s, 
unity. Without going further into the details of this process, which may be 
easily illustrated and worked out graphically, the results are as follows : — 

For n odd 



(T (e) = h s for values of e between 

' ^ n — \ y 2n J 
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and 
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(4) 



'■*For a given table n must not exceed twice the greatest tabular difference. 
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for values of e between and -I- ^ , (4 coni.) 



= — ■ e for values of s between H and + 

« — I I 2« J 2« 

For « even 
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= (i — e") for values of e between and i. 

n — \ ' n 

For n odd the probability area is a trapezoid symmetrical with respect to 
the axis of (f (e). For n even the area is bounded by three straight lines and 
the positive axes of s and (f (e), the errors as stated above being considered with 
respect to magnitude only. 

Now let Sy = the probable error, 

e„ = the mean error, 
e„ = the average error. 
The probable error is derived from the condition 

rfcSp 

I ^ (s) dfe = J for n odd, 
o 

= \ for « even. 

The square of the mean error 



sj ^2\ (f (s) ^ dt for n odd, 

o 
I 

= j ^ (e) sVs for « even. 



o 
The average error 

2« — ■ I 



2» 



E^-=-2\if(t) sds for n odd. 
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=j<f (£) ids 



for « even. 



In applying these formulae, attention must be paid to the several values of 
ip\{t) and the corresponding limits falling between those given in the above inte- 
grals. They give 

2« -f I I 



2« 



2« 



|/[2«(« — i)] 



(2« — l)' — I 

" 24«^(« — l) 

Sj, = I — ^ )/2 for « = 2 



> n odd ; 



(6) 



2« — 
2« 

for « > 2 



"~i'^[^^'^('^~')~^^ 



«« = 



«' + (« 



+ («-i)^-n 
■'(«-I) J 



> « even. 



2«'' 



i)^ 



6«^(« — i) 



(7) 



For obvious reasons these formulae, for n odd and even, should give identi- 
cal results for s^, e„, and e„, when « = oo ; moreover these results should be the 
same as those for n-=^2. Such is the case. The value « ^ oo is not permis- 
sible, however, in these formulae when applied to such tables as we use, since the 
tabular differences are always finite. To give consistent results when applied to 
such tables, n must be restricted to values not exceeding twice the greatest tabu- 
lar difference. For values of n greater than this limit, the probable, mean, and 
average errors will be less than those given by (6) and (7) and will approach 
those for the system expressed by equation (i) as n approaches infinity or unity. 
The process outlined above may be readily applied to determine the laws of 
frequency of s when n exeeeds the limit just mentioned. 

As to the more general case, in which the interpolating factor ^ is a proper 
fraction of the form w/«, m and n being integers, a moment's consideration will 
show that it is comprised practically, or so far as numerical values are concerned, 
in the first case where t = i/«. For, supposing m and n to have no common 
factor, the possible values of Sj corresponding to m/n are identical with those for 
i/n. Thus for m/n = f the values of f, are the .same as those for m/n = ^ ; 
those for ^ and -^ the same as those for ^ ; those for -f^ the same "as those for 



INTERPOLATED FROM TABLES OF LOGARITHMS, ETC. 



59 



yV, and so on. The values of f (s„) and hence (f (s) for / = mftt will indeed differ 
from those when t ^ ijn. <p (s) will be of a somewhat more complex character.* 
But the limiting values of s are the same in the two cases, and hence the 
probability functions cannot give probable, mean, or average errors differing ma- 
terially in numerical value from those when /^ i/«. The restriction mentioned 
above relative to the magnitude of n for a given table must be borne in mind. 

The following table gives the probable, mean, average, and maximum actual 
errors corresponding to the interpolating factor t = i/« for values of n from i to 
to. The values of these errors for n= i, or for a tabular value, are given in the 
first line for the sake of ready comparison. 
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The mean of the average errors in the above table for / = ^, ^, ^, .j^, 
and ^^, using -jlj and \ as arguments in the place of -^ and 3^, is 0.316. This 
number coincides closely, as it should, with the average of the lOOO actual errors 
given by Prof Howe in No. 6 of the Annals. 

*'rhus, for example, if / = Jj the probability function is 

(t) = I for values of e between o and ± -}^, 

-- J (ff ± e) for valBes of f between ± xV ^"d ± A> 
= j (I ± c) for values of t between ± -}^ and ± -}^, 
= J (xV ± *) fo"" values of e between dz j'j and ± t^ ; 
and the probable, mean, and average errors are 0.261, 0.351, and 0.290 respectively. 



